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Abstract. Using a new method of monotone iteration of a pair of smooth 
lower- and upper-solutions, the traveling wave solutions of the classical Lotka- 
Volterra system are shown to exist for a family of wave speeds. Such con- 
structed upper and lower solution pair enables us to derive the explicit value 
of the minimal (critical) wave speed as well as the asymptotic rates of the 
wave solutions at infinities. Furthermore, the traveling wave corresponding to 
each wave speed is unique modulo a translation of the origin. The stability 
of the traveling wave solutions with non-critical wave speed is also studied by 
spectral analysis of the linearized operator in exponentially weighted Banach 
spaces. 



We re- visit the classical Lotka-Volterra competition system: 



where u{x,t), v{x,t) are the population densities of two competing species, the 
constants di, ai, hi, Ci, i — 1,2 are assumed to be positive. In this paper, we are 
trying to acomplish the following goals: providing a new and easy construction of 
upper- and lower-solutions to derive the traveling wave solutions of (jl.ip ; obtaining 
an accurate description of the minimal wave speed and asymptotic behaviors (up 
to the first order) of the wave solutions; investigating the stability of the traveling 
wave solutions in various Banach spaces. 

System (|l.ip has been extensively studied. In |16-Leung] and [18-Pao| . there are 
many applications and treatments of solutions of (jl.ip in bounded spatial domain 
under various initial and boundary conditions. As is well known, system (jl.ip and 
its cooperative counter-parts admit traveling wave solutions. In |22-Tang| , Tang and 
Fife showed the existence of the traveling wave solutions connecting the extinction 
state with the co-existent state, In |09-Kane"I) , Kanel and Zhou studied the existence 
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of the traveling wave solutions connecting the coexistent state to a semi-exitinction 
state. In |05-Feij . the traveling wave solution connectiong two semi-extinction states 
were studied, they also estimated the minimal wave speed. For the other treatment 
of the traveling wave solutions of system (jl.ip and its generalizations, please see 
[25-Wuj . [01- Alexander] , [07-Hosonoj . |24-Volpert| . 



Throughout this paper, we make the following assumptions: 
. |H21. ^<^. 

C2 Ci 

• [H3]. + — < 1 + — . 

oi C2ai ai 

Under these conditions system (|l.ip has three equilibria (0, 0), (f^, 0) and (0, 

We will use a new monotone iteration method to investigate the traveling wave 
solutions of (jl.ip . The traveling wave solution connects two of the above equilibria. 
The monotone iteration method has been widely used in the study of the travel- 
ing wave solutions of reaction diffusion system such as (|l.ip , but most constructed 
upper- and lower-solutions in literature are non-smooth. The relatively larger 'gap' 
between the non-smooth upper and lower-solutions creates certain difficulties in 
deriving the accurate asymptotic estimates of the traveling wave solutions at infini- 
ties. Such estimates are valuable in applications, and enable one fully exploit the 
cooperative or competitive structure of the Lotka Volterra system. 

The smooth upper-solution in the monotone iteration as in section 3 is derived 
from the traveling wave solutions of the KPP equation. Observing that if we take 
fuction to be a constant, then the first equation of (jl.ip is a generalized KPP 
equation, the same consideration is also true for the second equation. The existence, 
uniqueness, asymptotics as well as the stability of the traveling wave solutions of 
the KPP system are well known, so are the properties of the upper-solution. The 
most difficult part in section 3 is to construct the lower solution for system (jl.ip . 
Since such constructed upper solution is 'nearly' a solution, a compromise is made 
to derive the smooth lower-solution. In fact, the lower-solution does not satisfy the 
boundary condition at oo. Thanks to the realxed condition in [25-Wuj . we can still 
apply the monotone iteration scheme as specified in |25-Wu| . [03-Boumenirj . The 
trade off of such 'shorter' lower solution is that we can have some freedom to choose 
the lower solution with desired asymptotic rate at negative infinity. This leads to 
an accurate (up to first order) asymptotic estimates of the traveling wave solution 
at — cx), and an exact value of the minimal wave speed. 

The asymptotics of the traveling wave solutions at infinities are obtained by 
comparison principle. Once the asymptotic behaviors of the traveling wave solu- 
tions are known, we can use the Maximum principle and Sliding domain method 
[15-Leung] to derive the uniqueness, strict monotonicity as well as the local stability 
of the traveling wave solutions. 

The local stability of the traveling wave solutions is studied by means of spectral 
analysis in some weighted Banach spaces. We proceed to show that the linearized 
operator about the traveling wave solution has essential spectrum lying completely 
in the left complex plane, and that is not an eigenvalue of the linearized operator 
in the weighted Banach spaces, whereas all the other eigenvalues of the linearized 
operator have negative real parts. This means the traveling wave solution is linearly 
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exponentially stable. Since such linear operator is sectorial, the linear stability im- 
plies the local nonlinear stability of the traveling wave solutions |01-Alexanderj . 
[24-Volpert] , |06-Henry] . Though general theories on the stability of the traveling 
wave solutions are known 24-Volpe"rt] , |12-Kapitula| , |01-Alexanderj . the verifica- 



tion of the conditions there is in fact a case by case study. The methods used in 
the stability study of the traveling wave solutions of are similar in the spirit 
to those in |28-Xuj . |21-Sattinger] , p-Kapitula] , [Ol-AlexandCT] , |27-Wu| . |26-Wuj . 
[02-Bates] , but in comparing to the above methods, we need to further overcome 
the 'unstable' component of the system. This is done by studying an equivalent 
form of the linearized operator in a smaller weighted Banach spaces. 

We remark that the stability is only for the traveling waves with non-critical 
wave speed. The stability of the traveling waves with critical wave speed is currently 
under investigation. 

The paper is arranged as follows: in section 2, we study the steady states of the 
system and obtain the attraction region of the stable steady state; in section 
3, we show there are traveling wave solutions connecting the one of the unstable 
steady states with a stable one, corresponding to each wave speed, the traveling 
wave solution is unique and strictly monotone. The analysis is done by utilizing 
an accurate description of asymptotic behavior of the traveling wave solutions. 
Furthermore, we have obtained the estimate of the critical wave speed. In the last 
section of the paper, we show the traveling wave solutions are locally, nonlinearly 
exponentially stable. 

2. THE EQUILIBRIA AND THEIR STABILITY 

In this section, we analyze the constant steady states of system (jl.ip under 
conditions II1-H3. 

Consider system with the initial conditions 

(2.1) u{x,0) = uo{x), v{x,0)^vo{x) xeR. 

where uq and vq are non-negative bounded smooth functions on R. 

Definition 2.1. A pair of nonnegative bounded smooth functions U = {u,v) and 
U = (u, v) are called coupled upper- and lower-solutions of the Cauchy problem 
p.ip - (|2.ip if [/ > [/ on Rx [0, oo) and the following inequalities are satisfied 



Ut 


> 


d\ Uxx 


+ u{ai 


— hiu 


- civ), 


Vt 


> 


d2Vxx 


+ v{a2 


- b2U 


- C2v); 


Ut 


< 


dlUxx 


+ u{ai 


— biu 


- civ), 


Vt 


< 


d2Vxx 


-\- v{a2 


- b2U 


~ C2V), 



(2.2) inRx(0,oo), 



as well as the initial conditions 

(2.3) m(x, 0) > uo{x) > u{x, 0), v{x, 0) > vo{x) > v{x, 0) in K. 

It is known from |16-Leung] , [18-Pao| that if there exist coupled upper and 
lower solutions U and [/ on R x [0,oo), then the Cauchy problem (|l.ip - (|2.ip has 
a unique solution U(x,t) — {u(x,t),v(x,t)) with u{x,t) > u(x,t) > u(x,t) and 
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v{x,t) > v{x,t) > v{x,t) on M X [0,oo). The next theorem gives the asymptotic 
stability and attraction region of the equihbrium (^,0). 

Theorem 2.2. Let a ^ hC^-T^). B ^ Mii(si±S2 _ al^ a ^= M2(ai_a2N 
Assuming that the hypotheses (HI) and (H3) hold, we know that a > Q, B > Q, 
and P > 0. If for some A > and < p(0) < ai/bi, the initial functions satisfy 

Q- - P(0), 0) < {u,{x), vo{x)) < + Ap{0),Bp{0)) 
on R, then the solution for (|l.ip - (|2.ip satisfies 

(2.4) - p(t), 0) < {u{t, x), v{t, x)) < + Ap{t),Bp{t)) 
on [0,oo) X R where 

(2.5) Pit) = [/3/a + (p(O)-i - /3/a)e"*]-'- 

Proof We will show that (m, u) = (f^ + Ap{t), Bp{t)) and (u, w) = (f^ - p{t),{)) are 
a pair of coupled upper-lower solutions defined in Definition (|2.1[) . One can easily 
see that -2 = satisfies the required inequalities in ()2.2|) . 

We first start with u — j^ + Ap{t). From (|2.2p . it needs to satisfy the differential 
inequality 

Ap'{t) > (^ + Ap)[ai -bii^+ Ap)] = -Aiaip + 6i V) 

Ol Ol 

For p{t) > 0, it suffices to show that 

(2.6) p'{t)+aip{t) > -biAp'^it). 

Also for V = Bp{t), we see from (|2.2p that it needs to satisfy the differential in- 
equality 

Bp'{t) > Bp[a2 - b^i^ -p)- C2Bp] 

Ol 

and it suffices to show that 

(2.7) P'W+^2(f^-ff)p(t) > {b2~C2B)p^{t). 

Finally we look at u = |^ — p(t). Again from l|2.2p . it needs to satisfy the differential 
inequality 

or equivalently, 



y(t)<(^-p)[ai-6i(^-p)-cii?p], 
Ol bi 

p'[t)>{^^p){ciB^bi)p. 

Ol 



For p{t) > it suffices to show that 

(2.8) p'{t) + f^{bi~ciB)p{t) > {bi-ciB)p\t). 
From hypotheses (Hf ) and (H3) we observe the fact that 

aib2 ,ai a2. 
ai > — 02 = 02(7 — ) = a > 

Ol Ol 02 

Setting f^(&i — ciB) = a leads to the choice of 

(2.9) B=^(^1±^_^)>0. 

aici 02 Ol 
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From the hypotheses (H3), one can obtain that 

(61 - ciB) ~ (62 - C2B) = (62 - — ) - (62 - C2B) 

ai 

61 C2 , , . b2C2 0261 

= (ai +02) 

aiCi Ci ai 

= — 1 + — - 7^ - ^-^ > 0. 

ci ai &1 aiC2 

Noting that bi — ciB — ^^(f^ ~ ff) = Z^: we can now conclude that the three 
differential inequalities (|2.6p - (|2.8[) will all be satisfied if the function p{t) is a positive 
solution of the differential equation 

(2.10) p'{t)+ap{t)=pp\t). 

This leads to the function p{t) given in (12. 5p with p(0) < a/j3 = ai/bi and 
linit^oo p{t) = . □ 



From the arbitrariness of constant A in Theorem (|2.2p . we then have the the 
attraction region for the equilibrium (ai/6i,0). When the hypotheses (HI) and 
(H3) hold, for all the initial density functions (uq^vq) in the rectangular area 

\ rr, ^2 + 0.2 fli 

0,00 X 0, — — — , 

Ci &2 Dl 

the solution {u,v) of the system (|l.ip - (|2.ip converges to the equilibrium (ai/6i,0) 
uniformly on R as i ^ oo with the rate e^"*. 

In the meantime, by adding hypothesis (H2), we can quickly find that the equi- 
libriums (0,02/02) and (0,0) are both unstable. For this purpose we construct a 
pair of upper-lower solutions 

{u, i) = (A/, a2/c2 - p{t)) and (m, v) ^ {Ap{t),0), 

where M > ai/bi is a constant. Constant A and function p{t) will be determined 
later. The differential inequalities in (|2.2p are automatically satisfied by u and v. 
For u and v, the following relations need to hold: 

p'^t)<iai~^)p+ic,~Ab,)p^ 

C2 

and 

P'(i)<(--P)(c2-A62)p. 
C2 

The above inequalities are equivalent to 

p'it)-ia,-^)p<-{b^A-c^)p^ 

(2.11) 

p'{t) ~ ^{b2A- C2)p < -{b2A- C2)p^. 



a2Cl 0,2, , 

Oi = (62A - C2), 

C2 C2 



Setting 

from hypothesis (H2) we have 

(2.12) A=l(c2-ci + ^)>f, 

02 02 02 
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From the hypotheses (H3), one can obtain that 

[biA - ci) - (62^ - C2) = 7-(c2 - ci H ) 

62 0.2 0,2 

= (ai + 02) 

a\C\ c\ a\ 

_ ai6iC2 j-j^ _|_ f^. _ ^ _ «2Ci j ^ Q 
0262 ai &i aiC2 ~ 
Both the incquahtics in (|2.1ip can be satisfied by choosing p(t) as the solution 
of the differential equation 

(2.13) - 7P = -<^p' 

where 

7 = fli ^— !- > and 5 = h\A — c\> 62^ — C2 = ^ ^ — ci > 0. 

C2 02 

This results in the function 

with an arbitrary constant C > 0. For arbitrarily small e > 0, one can always 
find a constsnt C large enough such that p(0) = 7/((5 + C) < e. The fact that 
limt^oo = 7/'^ leads to the following theorem indicating that (0, a2/c2) and 
(0, 0) are both unstable. 

Theorem 2.3. Let 

(2.15) 7 = ai ,A^—{c1-c\^ ) and d = — (c2-ciH )~c\. 

C2 02 02 02 02 

Assuming that the hypotheses (HI), (H2) and (H3) hold, we know that 7 > 0, 
A > C2/62 and 5 > — Ci > 0. For any arbitrarily small e with < e < 
min{747/(5, 7/(5}, if the initial Junctions (uo,fo) satisfies Uo{x) > e and < vo{x) < 
ailc-2 — e on f/ien i/ie solution {u{x,t),v{x,t)) of (|l.ip anrf (|2.ip satisfies 

A^ fl2 7 

(2.16) liminf u(a;, i) > — ^ and limsupu(x,t) < r. 

t^oo t-»oo C2 

3. THE TRAVELING WAVES 

In section [21 we showed system (jl.ip has two unstable constant steady states: 
(0,0), (0, ||) and one asymptotically stable constant steady state (|^,0). We will 
show that there are traveling wave solutions of (jl.ip having the form 



(3.1) {u{x, t), v{x, t)) = {kw{\j -jx + cait), qz{\^ -jx + coit)), 

and connecting the unstable state (0, ||) with (^7,0) as the variable y^x + coit 



runs from —00 to +00. The constant c in (|3.2p is the wave speed and the minimal 
speed is also called the critical wave speed. Throughout the rest of the paper, we 
assume di = ^2 = d. 

To simplify notions, we introduce the following transformations to 
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(3.2) 6 = 62fe^^ £2 = a2^^C2g - 1, 

k — aibi^ , and 5 is a constant satisfying a2C^^ < 9 < aic^^. 
Under transformations (|3.ip and (|3.2p . system (jl.ip is changed into 

{— + cifj = w{l — w — rz) 
-cx)< ^ < 00, 
— + cz^ = z(ei — 6w — ei(l + 62)2) 

with the corresponding boundary conditions 



(3.4) 



1+62 

= (1,0). 



where ^ = V?'^ + '^"1* ™ <E31), (E3) for a; G M and t e 
We further introduce the transformations 



(3.5) uiiO^wiO, U2{0 = T^^<0 

to change system p.3p into the foUowing monotone (cooperative) system: 



(3.6) 



f \ / N A + <^2 — r 
-(ui)« + c(ui)5 = ui[—— ui+ru2), 

i + £2 

-(""2)5^ + 0(7/2)4 = (— -^^ U2)(bui - ei{l + e2)u2) 

1+62 



with boundary conditions 

hm4^_oo(«;(e),^(C)) = (0,0), 



(3.7) 



1+62 



Remark 3.1. Note that from hypotheses [H1]-[H3] and relations (3.2), we have 
the folfowing inequahties: 

< ei < 5, < r < 1, £2 > and 1 - > 6 - ei > 0. 

Before showing the existence of the travehng wave solutions for (|3.6p with bound- 
ary conditions (|3.7p . we first recall the following well known fact: (please see 
14-Kolmogorov] , [21-Sattinger| for the proof) 



Let a function / be a function on the interval [0, /?],/? > 0, with / > on 
(0,/3), and /(O) = /(/3) = 0, /' (0) = ai > 0, /'(/?) = -/3i < 0. 
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Lemma 3.1. Corresponding to every c > 2,^fa\, the boundary value problem 

(3.8) { e e M 

tt'(— oo) = 0, uj{+oo) — b. 

has a unique monotonically increasing traveling wave solution ti-'c(C) j C S M , where 
the lower index denotes the dependence of the wave solution to on c. 



We next show the existence of the travehng wave solution for system (I3.6l) - (|3.7p . 

Theorem 3.2. Let the parameters ei, £2, b and r satisfy conditions in Remark 3.1, 
then corresponding to every c > 2^1 — system US. 6\) has a monotone traveling 

wave solution satisfying the boundary condition |j'.7| ). (Recall hypotheses [H1]-[H3] 
imply all the conditions in Remark 3.1 are valid.) 

Proof. The proof will be done by monotone iterating a pair of smooth upper- and 
lower-solutions. We first construct a twice differentiable smooth upper-solutions. 



According to lemma [331 foi' every c > 2^1 — y^t^i there is correspondingly a 
montonically increasing function F(<^), C G K satisfying 



(3.9) 



r(-c») = o, y(oo) = 1. 



Define 



(3.10) MO=y{0 , U2{i)^-^Y{i). 

1+62 

For < U2{£) < U2{£), we readily verify that 

-u'iiC) + cui(C) - uiC-^ - ui -f ru2) 

I + £2 

r _ _ _ r _ 
= (1 - -r- - wi) - ui(l - — ui+ru2) 

1 + 62 1 + 62 

(3.11) = "i{(i-TT^)(i-"i)-i + TT^ + "i-™2| 

_ r r _ 

= ui < — ui - ru2 

L 1 + £2 

_ , r _ 
> ui { ui - ru2 

1 + 62 
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for all — oo < ^ < oo. For < ui{£_) < ui(^), one verifies 

-u'^iO + cu'^iO - (t^ U2){bui - ei(l + £2)^2) 

^ {-Y" + cY' - {1 - Y){bui ~ eiY)} 



1 + 62 



(3.12) = + 

^ -(i-y)r(i--^ + ei-4 >o. 



1 + 62 I 1 + 62 



The last inequality is true provided b < 1 — y^t^+Ci, which is valid due to hypothesis 
H. It is also straightforward to verify that (ui, U2) satisfies the boundary conditions 

We next construct a twice continuously differentiable lower solution for the system 
dS^-dnZl). Let the function Z(^), ^ e R be the solution of 



(3.13) 



% - cZ^ + (1 - ^)Z(1 - T3zg^Z) = 0, 

1 + ^2 



Z(-oo) = 0, Z{oo) = J^^^ 



Here I is some number in the interval (0, 1) to be determined. One can readily 
verify that the solutions of (|3.9p and p.l3p are related by the following 

(3.14) z{0^ ]l^ Y{0, eeK- 

1 + 62 

Since < ^ < 1, we have 

(3.15) Z{0<Y{0, eeM. 
We define a lower solution of (|3.6p . (|3.7p by setting 

(3.16) ui = Z, U2 = Z. 

1+62 

where I G (0, 1) is to be determined. We readily verify that they satisfy 
-"1 (0 + cui(0 - Mi( t ,^ ui+ rU2) 

(3-17) ^ z((l-^)-(l-^)Z-(i±i^) + Z-^4 

V 1+62^ ^ 1 + 62^ ^ 1+62 ' 1 + 62 J 



= 0. 
Moreover, we have 
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(3.18) 



+ cu'^iO - i-T, U2)ibui - ei(l + 62)722) 

1 + 62 

Z i 1 1 )Z } Z) {hZ - eilZ} 

1 + 62 I 1+62 ^ l + e2' j 4+62 1 + £2 ^ ^ ^ 

^ -z\{l-Z)- r(-i - (— T-^^) - ei^^} 



1 + 62 I 1 + 62 1+62 J 1 + 62 1+62 

< (- Z)\-^—Z-bZ + eilz\ + il-lZ)- 



1 + £2 1 + £2 I 1 + £2 J 1 + £2 

= —Z)l —-b + e,l + l\z 

4 + 62 l+£2 l+£2 J 

< 0. 

The last inequality is valid provided that 

(3.19) 0<Z<niin{l,-; —}. 

Also by the limiting boundary conditions of (|3.13p we see (mi.U2)(— oo) = (0,0) 
and {ui,U2){+(x,) = (ttCT' (T^7)(TTi^)- Inequality along with (EH]) 

show that (wi,M2) consists of a pair of lower-solutions for system (|3.6[) . (|3.7p . 

Noting that such constructed upper- and lower-solution pairs are ordered. We 
can apply the monotone iteration methods provided in |25-Wu| or |03-Boumenir| 
to derive the conclusion of this Theorem. Here we only sketch the ideas. 

Following the notions in |25-Wu| . we write (3 = diag. (0,0), K = {Ki,K2) := 
(li iq^). lower-solution p(^) = (Z(^), j^Z{^)), and upper-solution p(^) = (r(^), 
j^^y{£,))j £ as described above. As explained below, a slight variant of The- 
orem 3.6 or Theorem 3.6' in [25-Wuj is needed, because {Fi,F2) defined in p.22p 
below has an additional zero at (0, j^) between := (0,0) and K. By means of 
the iterative procedure in the proof of Theorem 3.6 in |25-Wu] . we first obtain a 
solution of (|3.6p (j){£_) :— (ui(^), U2(C))i satisfying the inequality 

(3.20) m<m<pio, ^eK- 

From the comparison argument with p{^) in the proof of Theorem 3.6 in [25-Wuj . 



we 



have 



lim Mi(^) — lim W2(0 = 0. 

^ — * — oo ^ — > — oo 



Again, by the limit and comparison argument in the proof of Theorem 3.6 in 
|25-Wu| ■ we obtain for i = 1, 2, 

(3.21) lim u,(0 = Q^, F,{Qu Q2) = 0, 

where 

p, I r., — r,, (I 

1 + C2 



Fi{pi,P2)=Pim^-pi+rp2). 



(3.22) 

^2(Pl,P2) = (l^ - P2){hpi -£1(1 + £2)^2); 
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and 
(3.23) 

< lim Z(e) < Qi < i^i - 1, < lim -J— Z(e) < Q2 < K2 = 

C^oo ^^oo 1+62 i + £2 

We then deduce from condition (H) and (I3.20p . p.2ip . (|3.23p that we must have 
(3.24) Qi=i^i = l, Q2 = K2 ^ 



1 + e2 

□ 

Remark 3.2. One can translate Theorem 3.2 into the foUowing: Assuming hy- 
potheses [HI], [H2], [H3], then for any c > 2^1 - system pTT|) has a travehng 

wave solution connecting (0, ^) with (|^, 0) as the variable ^/^x + cai< running 
from — CX3 to -\-oo. 

Remark 3.3. Theorem 3.2 does not insure the strict monotonicity of the resulting 
traveling wave solutions, as the iteration is only applied for the upper-solution, and 
the lower solution is served as a nonzero barrier so that the iteration limit does not 
coverge to zero. 

To further study the asymptotics of the traveling wave solutions as obtained in 
Theorem 13.21 we shall need the following Lemma concerning the scalar problem 
(3.8). 

Lemma 3.3. The solution Wc{S.) to (3.8), described in Lemma 3.1, has the following 
asymptotic behaviors: 

1. Corresponding to the wave speed c > 1JlSi\, 



(3.25) 



^c(c) = a^e 2 « -I- o(e 2 as 4 ^ -00 

^c(C) = P - Oi^e 2 '^+o{e 2 as +00, 



where a^j and b^j are positive constants; 

2. Corresponding to minimal wave speed c = 



(3.26) 



Wc(0 = (flc + 4e)e^« + o{S,e^^), as <e -00 
iOc(^) ^ P - bee 2 «-|-o(e 2 t), as ^ 



where the constant dc < 0, be > and Oc G M. 

Proof. The conclusion follows |21-Sattinger| , |23-Thiery] with slight changes. □ 



Based on Lemma 13.31 we study the asymptotic behaviors of the traveling wave 
solutions of system (13.61) . p.7p at infinities. 

Corollary 3.4. Assume hypotheses [HI] to [H3], and thus all the conditions in 
Remark 3.1 are satisfied. Let a = 1 — yipjj, then the traveling wave solutions 
(ui(^), M2('C)) of system i fg.6)) . \3.T^ as obtained in Theorem \3.2\ have the following 
asymptotic behaviors: 
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1. Corresponding to each wave speed c > the traveling wave solution 

(ui(^), 1*2(0) satisfies 

(3.27) U + o(e^^«) 



"2(0 / \A 



2 



as ^ ^00; and 



MO \ f ^ \ f 



(3.28) 1 1 = 1 ^ 1-1 |e ---^^^t^ ^+o(e ---^^^?^ ^^ 

.2(0 ; [^j ' 

as ^ ^ +00, where Ai, A2, Ai, A2 are positive constants; 

2. Corresponding to the wave speed Cqj-iUqqI = "^.^/a, the traveling wave solution 
(ui(0,'"2(0) satisfies 

( Ul(0 \ / All, + Al2rX \ 

(3.29) = ev^« + o(^e^«) 

V U2(0 / V ^21c + A22c^ ) 

as ^ —I- —00; and 

,3.30, I "'IJI ) = ( ^ ) - ( J" ) .-^^=.„,.-^^=) 

as ^ ^ +00, where Ai2c, ^22c < 0, Anc, j421c G R anrf ^11,^22 > 0. 



Proof. For c = 2y^l — jipj;-: according to Lemma 1331 the upper-solution (ui,U2) 

and the lower-solution ('Ui,-U2) as defined in (|3.10p . (|3.16p have the following re- 
spective asymptotic behaviors at —00, 



WlO \ / Aiic + Ai2c^ 

U2{0 ) \ A21c+A22c^ 

and 

^^2(0 / V B21c+B22c^ 



3V^? + o(Ce^«), 



While with wave speed c > 2^1 — 

^^2(0 / \A2 



c-Vc2 -4c, ^ c-Vc2-4c. ^ 

e 2 « -I- o(e 2 «)^ 



and 
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/ \ ( M 

V "2(0 / V ^ 

where Ai, Ai, ~Ai2c, — i?j2c 1 — 1,2 are positive constants and Anc, Buc are 
constants. 

Noting the upper- and lower-solutions have the same asymptotic growth rate at 
— oo, we immediately have p.27p and p.29p via comparison. 

We next derive the asymptotic behaviors of the traveling wave solutions of (|3.6p . 
p.7p at -l-oo. Since the the traveling wave solution {ui{(,),U2{^))'^ monotonically 
aproaches the steady state (1, jqij^) as ^ — > -|-oo, by letting (wi(0,W2(0) be the 
derivative of the traveling wave solution, we have (wi, W2)(4-oo) = (0, 0) |25-Wuj . 

We linearize system (I3.6P about the traveling wave solution {ui{^),U2{£,)) to 
obtain 



c- Vc^ -4 



o(e" 



{055 - + (1 - if^ - 2iii + ru2)<l) + ruiip = 0, 
V-CC - CT/^? + b{j^^ - U2)(f> + [-ei - bui + 2ei(l + 62)7/2]^ = 0. 

It is easy to see that {wi, W2) solves p.3ip . 

The limit system of p.Sip at +00 is a constant coefficient system, and is given 

by 



(3.32) 



(0+)?e-c(0+)^-0+ + rV^+ =0, 
(^+)??-c(^+)? + (ei-W+ = 0. 



The exponential growth rates of the traveling wave solutions of (|3.6p - (|3.7p are 
determined by those of the solutions of p.32p . The justification is as follows: first 
we note that (|3.32p admits exponential dichotomy. By the roughness of exponential 
dichotomy, solutions of p.3ip grow/decay exponentially (possibly with a different 
exponential rate) |04-Coddington| , [20-SandstcdeJ. Since the derivative {wi,W2) 
of the traveling wave solution solves ()3.3ip . the traveling wave solutions of ()3.6p 
approach exponentially to the steady state (1, j^^)- 

To find out the exact asymptotic rates of the solutions of (|3.3ip , we first change 
p.3ip . (|3.32p into the first order systems, letting cj)^ = 4'^ = "ipi, <l>^ 
, we have 



-I- 
1 ' 



d 



/ \ 




( ^ \ 


01 




01 






\ ^1 ) 




^ V'l / 



and 





( 0+ 


] 




/ 0+ 




d 












dl 














) 






) 
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where 



7^^ 



and 





-l + ^ + 2u,{0-ru,{0 c 







1 \ 

-rUliO 

1 

ei+6wi(0-2ei(l + 62)1x2(0 c / 



/ 1 ^ 

1 c -r 

1 

yO 6-ei c/ 

The exponential growth of the travehng wave solution (ui(0,U2(0) at infinity 
implies that 

POO 

/ |7^+ - 7^°° < +00 for some large to > 0. 

J to 

It is easy to check that the matrix TZ°° has 4 distinct eigenvalues, then by [04-Coddington| and 
|21-Sattinger| , the asymptotic exponential rates of the solutions of p.3ip are the 
same as those of the solutions of p.32p . 

We search for the solutions of (|3.32p with zero limit at +00. Owing to conditions 
in Remark 3.1, the second equation of (|3.32p has general solution of the form: 



(3.33) 



-Vc^+4(6-8i) , 



C22e 



c + Vc^+4(b-ti) , 



with 



<o, 



> and C21, C22 two constants. 



Substituting (|3.33[) into the first equation of (|3.32D . we then have 



(3.34) (0+)«-c(0+)«-0+ = -r(c2ie- 
with general solution of the form 



-Vo^+4(i,-ei) , 



C22e 



c+Vc^+4(b~-ti) , 



(3.35) 



ci2e 



C2ie 



-Vc'-'+4(b-ei) , 



C22e 



e + Vc^+4(ti-»i) , 



Comparing (3.31), (3.32) and using the fact that (wi, W2) is a solution of (3.31) 
with (wi(+oo), W2(+oo)) = (0,0), we deduce from (3.35) that as ^ ^ +00, 



wiiO = (ci2+o(l))e' 



(521 +o(l))e- 



-Vc^+4(t.^ei) , 



W2iO = (C2i +o(l))e 



= -Ve^+4(6-i!i) , 



with C21 ^ and 621 ^ 0. Since ^ 
write 



< 0, by [H3], we may 



c-Ve^+4(ti-ei) , 



wi(0 = (521 +o(l))e- 
as ^ ^ +ao. 

We therefore conclude this proof by integrating (wi(0, W2(C))"^ from ^ to +00 
with ^ sufficiently large. □ 
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We next derive a corollary that is important to section [H 
Corollary 3.5. Assume the hypotheses of Corollary 3.4- For every wave 
with c > 2y^l — j^i^, the corresponding traveling wave solution (wi(^), M2(^))"^ of 
system 13. 6\} . Tj j obtained in Corollary 3.4 is a strict monotonic function for 

C e R. 

Proof. According to the monotone iteration ( |25-Wu| ). the traveling wave solution 
('Ui(^), ■U2(C))"^ to (13. 6|) . (13. 7p is monotone, which implies that ""2(0)"^ ^ 

for ^ e R. By Corollary [331 one sees that «(±cx.), u'2(±oo))^ = 0. The 
monotonicity of system (j3.6p and the maximum principle lead to the conclusion 
that K(e),'a'2(C))^ >0 for CeM. □ 

We next show the uniqueness of the traveling wave solutions. 
Theorem 3.6. Assume hypotheses [HI] to [H3]. The traveling wave solution to 
system i3. 6]}-[37?\), obtained for each wave speed c > 2^ /I — 



-, with properties 

described in Corollary 3.4 and 3.5, is unique up to a translation of the origin. 

Proof. We only prove the conclusion for traveling wave solutions with asymptotic 
behaviors (|3.27p and (|3.28p . since other case can be proved similarly. Let C/i(C) 
and U2{£,) be traveling wave solutions of system ()3. 6^ - 1)3. 7p . with speed c > 
and properpties described. There exist positive constants Ai, Bi, i = 1,2,3,4 and 
a large number > such that for f < — A, 



(3.36) 
(3.37) 

and for ^ > N, 
(3.38) [/i(0 



(^1 +o(l))e- 



V {A2 + o{l))e' 



(A3 + o(l))e ^'^T^ C 



{Ai + o{l))e ^''^f^ i 



l-(Bi+o(l))e- 



-Vc^+4(b-ei) 



1 , , c-Vc^+4(i,-.i) 

(B2 + o(l))e 2 ? 



1+62 



/ l-(B3 + o(l))e 



:-Vc^+4(b-ti) 



(3.39) 



U2iO = 



- (Bi + o(l))e 5 ? 



V 1+62 

The traveling wave solutions of system (j3.6p are translation invariant, thus for any 
9 > 0, U^iO ■= Ui{^ + 9) is also a traveling wave solution of ([SJ]) . By ([3^ and 
()3.38p . the solution Ui{£, + 9) has the following asymptotic behaviors: 

/ 

(3.40) UliO = 



(^i+o(l))e 



:+Vc^-4q a c+Vc^-4q 

2 "e 2 



V (^2+0(1)) 



c + Vc2^4q a c+Vc2-4q 

e 2 "e 2 
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for ^<-N-9; 
(3.41) U'.iO 



{Bi + o(l))e 5 f>e ^ — 



V IT 



£2 



- (^2 + o(l))e 2 



for e > TV. 

It is clear that for 9 large enough, we have 



Aie' 



A2e 



-Vc3+4(b-ei) , 



> ^3, 



^26 ^ ^ < B4. 



(3.42) 
(3.43) 
(3.44) 
(3.45) 

Thus for some iV > 0, formulas (|3.42p - p.45p imply that for 6 large enough, 
(3.46) C/f (0 > C/2(C) 

for ^ e (— cx), — iV]U[Af + 00). Here, the inequality " > " in (|3.46p is component- 
wise. We now consider system (13.61) on [— A^, +iV]. First, suppose Uf{^) > U2{C) 
on [-N,+N], then the function W{Cl := U^{Cl - 1/2(0 > and satisfies for some 
CG (0,1): 
(3.47) 

f^i(C/2+Ci(t/f-C/2)), it(C/2 + Ci(^f-f/2)) 



W" - cW + 



if (C/2+C2(t/f-C/2)), ft(C/2+C2(t/f-C/2)) 



dF2 , 



W = 0, 



for ^ e (-iV,iV), and W^(-iV) > 0, W{+N) > 0. 

Since the above system is monotone, we can readily deduce by maximum principle 
that > 0, on [-N,N]. Consequently, we have J/f (C) > C^2(C) on R in this case. 

Second, suppose there are some points in (— TV, iV) such that [/f(0 < C^2(C)- 
We then increase 9, that is shift /7f (^ further left so that U^{-N) > U2{-~N), 
Uf{N) > U2{N). By the monotonicity of C/f and 1/2, we can find a 6* e (0,27V) 
such that in the interval (-iV,iV), we have Uf{£, + 9) > U2{£.)- Shifting Uf{£, + 9) 
back until one component of J/f (^ + 0) first touches one component of [/2(C) at some 
point ^ € (—-/V, N). Then by maximum principle for that component again, we find 
that component of Ui and U2 are identically equal for all ^ G [—N^ N] for a larger 
9 than the original one such that (|3.46p holds. This is a contradiction. Therefore, 
we must have 

uliO > U2{0 

for all ^ E R, where 9 is the one chosen by means of p.42p - p.45p as described 
above. 

Now, decrease 9 until one of the following situations happens. 
1. There exists ^ > 0, such that C/f (^ = t^2(C)- this case we have finished 
our proof. 
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2. For ^ > 0, there exists G R, such that one of the components of and 
U2 are equal at the point ^i; and for aU ^ e M, we have Uf{£,) > C/2(0- We then 
consider the system p.47p on {—N, N) and 6* = ^ in the definition for W. To fix 
ideas, we suppose that the first component of Uf and U2 is equal at the point ^1. 
The maximum principle for this component implies that the first component of 
C/f (^) is identically equal to that of U2{0- Also, we readily obtain that for large 
the limiting equation for (|3.47p is the same as p.32p . Since the first component 
of W is identically zero and the off diagonal limit coefficient r in the first equation 
in (3.32) is not equal to zero, we conclude from the first equation in (|3.47p that the 
second component of W must vanish for all large ^. By the maximum principle for 
the second equation, we conclude that the second component of W is also identically 
zero for all ^ G R. Similarly, we next consider the case that the second component 
of Ui and U2 are equal at the point ^1. We first obtain the limiting equation of 
(3.47) for large — ^. The off diagonal limit coefficient of the second equation will be 
^ 0. We first deduce that the second componet of W is indentically zero, and 
then the first component must also be identically zero for ^ G K. 

Consequently, in either situation, there exists a ^ > 0, such that 

ukO ^ U2{0- 

for aU C e M. □ 

Theorem 3.7. Assume hypotheses [HI] to [H3]. System p.6[) - (|3.7p does not have 
strict monotonic traveling wave solution tending to (0, 0)"^ as ^ — > —00 for c < 2y/a. 
Here, a = 1 - jf^. 

Proof. Suppose there is a constant c with < c < 2y/a and a corresponding solution 
V{^) = iviiO,V2iOV of (EH) tending to (0,0)^ as ^ ^ 00. Similar to the proof 
of Corollarv l3.4[ we can deduce by integrating the asymptotic approximation of its 
derivative that the asymptotic behaviors of V{(,) at —00 must be of the form: 

MO \ f As\ ( A 



e 5 5 +1 e 2 5 + h.o.t, 

V2iO J \Bs J \Bs J 

where {Ag, Bs)^and {As,Bs) can not be both zero, and h.o.t. is the short notation 
for the higher order terms. The condition < c < implies that V{£,) is 

oscillatting. This says that such solution of (|3.6p with c < 2^/a is not monotone. 

□ 

4. STABILITY OF THE TRAVELING WAVES WITH 
NON-CRITICAL SPEEDS 

In this section, we always hypotheses [HI] to [H3] for system (1.1); thus all 
the conditions in Remark 3.1 are satisfied for (3.6) and subsequent systems. We 
first show that the traveling wave solutions with the non-critical speed obtained in 
Theorem l3.2l is unstable in the space of continuous function C(R) x C(R) (Please see 
definitions below). This motivates us to investigate the stability in the "smaller" 
exponentially weighted Banach spaces. We will concentrate on the stability of the 
traveling waves with non-critical wave speeds. 

The following set up of the problem is standard: Let U = {u, v)^ , F{U) = 
(^i( ^t+£2'' ~ "1 + ™2), ij^ - U2)[bui - ei(l + £2)^2))^ and write system ([5^ 
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in the moving coordinates £, — x + ct, c > 0. In terms of t) variable, ()3.6p is 
changed into 



(4.1) 



f/(e,o) = U, 



where U the initial function. Let J7*(^) = {u* {^) , v* {^))'^ , ^ = x + c*t he the 
traveling wave solution of (|3.6p with wave speed c* > l^fa . It is easy to see 
that t/*(^) is a non-constant steady state of the system (|4.1[) . We consider the 
perturbation of this steady state solution. By letting U{£,,t) = U*{£,) + V{^,t) , we 
obtain the system for the perturbation function V: 



Vt = LV+U{V,U*), 



(4.2) 
where 

(4.3) LV = V^^-c*V^ + ^{U*)V 
is a linear operator, and 

(4.4) JV{V, U*) = FiU* + V)~ FiU*) ~ ^{U*)V 

is a nonlinear operator. 

The stability of the traveling wave solution C/*(^) in certain Banach space is 
determined by the location of the spectrum, (t(L), of L. 

Let 

o'p{L) = {A e cr(i) \X isan eigenvalue of L}, 

and ae{L) be the essential spectrum of L, which are points in a{L) outside <Jp{L) n 
{isolated eigenvalues of L with finite multiplicity} . Note that ae{L) includes 
the continuous spectrum of L. |01-Alexanderj . |06- Henry] , |24-Volpert] . Let C(M) 
be the space of all continuous functions on the real line and Co(K) be its subspace 



Co(M) := {C/(0 e C(R) X C(M) I lim U{£,) = 0} 



along with norm 



U llco(R)— sup II U{0 
5eR 



We also need the following weighted Banach spaces: for non-negative numbers cti, 
(72, the space Ca-i.a2 is defined as: 

= {U{0 e Com I UiOie'^'^ + e-"^«) G Co(M)}, 
on which we define the norm 

II lie, .,= sup|| (7(e)(e"^« + e-"^«) || • 

Similarly, we can define Caua2: * = 1; 2, ... as well, for example: 

Ci%2 = {U\ UiO, U'iO, U"iO e C.,,.^, and C e M} 
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with 



m 



sup S,^o 



It can be readily verified that these spaces are Banach spaces. 

Theorem 4.1. Assume [HI] to [H3], and let U*{^) = [uiiO^MOV be the trav- 
eling wave solution of (3.6) or (4-1) with wave speed c* > 2^/a, a = 1 — j^:^ as 
described in Corollary 3.4- Then unstable with initial conditions in Co- 



Proof. This theorem holds for the traveling wave solutions with critical and non- 
critical wave speeds. We need to prove the trivial solution of (14. 2|) is unstable. 
Thus, it suffices to show that in the space Cq the operator L in ()4.3|) has essential 
spectrum with positive real part. As is well known ( [06- Henry] , [24-Volpert] ) the 
location of the continuous spectrum of the operator L is bounded by the spectrum 
of L at ±oo, which we denote by and L~ respectively. More precisely, we let 



L+V 



(4.5) 



V^i-c*V^ 



-1 

ei 



V, 



L-V 



(4.6) 



1 - 



1 + 62 



b 





±oo. 



V. 



Here, Uj. respectively denote the limit of U*{^) as ^ 
Now consider the equation 

Following |24-Volpert] and [06-Henry] , we replace V by e'^'''*+*''^)/, where / is an 
identity matrix and A is a complex number and C is real. We then have 



(4.7) 



{^C'i~c*Cti + —{ui)~xi)^o. 

The spectrum of the operator consists of curves given by: 



(4.8) 



BF 

det{~ei - c*CtI + —([/;) - A/) = 



Solving (|4.8p . we have 

(4.9) 

or 



dU 



C - cCi - 1 - A = 0, 



(4.10) -C^-c*C« + ei-&-A = 0. 

Letting X — x + yi ior x , y E M., then by (|4.9p we have 
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or by gUDl), 

(4.12) ^ = + 

Similarly, the spectrum of consists of curves: 

or 

(4.14) a: = -^-61 

in the complex plane. Consequently, by theory described in [06-Henry| , we have 

r r 
maxi?e(Te(L) > max{— 1, ei — 6, 1 — , — ei} = 1 — -; > 0. 

Hence, by |06-Henry] again, the traveling wave solution U*{£,) of (|4.ip is essen- 
tially unstable in Co(K). □ 

In order to obtain stability for the traveling solution U* , we will restrict the initial 
conditions and the operator i to a "smaller" Banach space Co-i,ct2 with cri > 0, 
(72 > and af + a2 ^ 0. To relate the operator L in Co(K.) to an equilvalent 
operator in CrTi,a2i we introduce the mapping T : Caia2 ^ C'o as follows: 

(4.15) TV := {e"^^ + e'''^^)V. 

T is thus linear, bounded and has a bounded inverse : Co — > Cai,a-2 with 
T^iy = (e'^i« + e"'^2?)-iy. Consider operator 

(4.16) L:V ^ TLT-^V. 

One readily sees that L is a linear operator with domain C'^^^{M) x C'^^^M). By 
relation ()4.16|) . considering L in Co-j^o-j is equivalent to considering L in Co(M), 
which is: 

dF 

(4.17) LV = - (251 + c*)V^ + {2gl -32 + c* + q^{U*))V, 

where {2gi+c*), 2gf—g2+c*gi in (|4.17p are short notions for the matrices {2gi+c*)I 
and M(^) ^ {2gf — 32 + c*.gi)/, where 



51(e) 
52(0 



fi"^!? -I- R-°'2? ■ 



with 



lim^^oo 51(C) = 0-1- lim«^-oo5i(0 = ^^ra; 



lim^^oo 52(C) = lim^^_oo 52(0 = cr|. 

We now locate the essential spectrum of the operator L in the space Co (I 
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Lemma 4.2. Suppose oi and ai satisfying 



(4.18) 



< (Tl < ^ 2 ' 

< ^ 5^ iii^ < (72 < ^ T^— 



then the essential spectrum of the operator L in the space Cq(R) is contained in 
some closed sector in the left half complex plane with vertex on the horizontal axis 
left of the origin. Outside this sector, there are only a finite number of eigenvalues 
ofL. 

Proof. As in the proof of Theorem 14. 1[ we first study the opearator L at infinity. 
We have 

(4.19) L+V = y« - (2ai + c*)V^ + (a? + cVi + -g^{U*))V, 



(4.20) 



OF 

L-V = % - (-2(72 + c*)V^ + (4 - 0*02 + Q^iU*))V, 



where a{ + c*ai + ^{U*) and — c*(72 + ^(U*) correspond respectively to the 
matrices: 



(4.21) 

and 

(4.22) 



M+ = 



M- = 



rf +cVi - 1 


- c*a2 + 1 



al + c*(7i + ei 



1 + 62 



b 

l+e2 



C (72 - ei 



Similar to the proof of Theorem 14.11 we find the right most points of the corre- 
sponding parabolas are on the horizontal axis given by 



(4.23) max .{(7?+c*(7i — 1, af+c*ai+ei—b, cr^—c* (72+1 , erf — c*(79 — ei}. 

1 + 62 

A simple calculation shows the number above is negative by the choice of cri and 
(72 in (4.18). Thus by the theory in |24-Volpert] , the essential spectrum of L is 
contained in a closed sector in the left complex plane with vertex on the horizontal 
axis left of the origin. Moreover, we may choose this sector with the further property 
that outside it there is a finite number of eigenvalues of L. □ 

Corollary 4.3. Assuming all the hypotheses of Lemma \4-.S\ and ai, satisfying 
\^.18^ , the essential spectrum of the operator L in the space Ca-i.a-2 contained in 
some closed sector in the left half complex plane with vertex on the horizontal axis 
left of the origin. Outside this sector, there are only a finite number of eigenvalues 
ofL. 



Proof. The conclusion follows immediately from Lemma [4.2l and relation (|4.16p . □ 
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Having established the location of the essential spectrum of the operator L in 
the space Co-i .0-2 j we next study the location of its eigenvalues. We first note that 
from Corollary 3.4, for c > 2^, we have {U* {£,))' {e'^^^ + e"'^^^) is unbounded as 
^ ^ — 00, which is different from the situations met in [02-Bates] . [28-Xu| . therefore 
their methods can not be carried over to our case. 

Lemma 4.4. Let ui and (J2 satisfy (4-18). Then is not an eigenvalue of the 
operator L in the space Ca-ia2(M)- 

Proof. Let U* be a traveling wave solution of (|4.1[) as obtained in Theorem 13.21 It 
is easy to see that (U*)' £ Cq and satisfies the equation 

(4.24) LV = 0. 

This shows that is an eigenvalue of the operator L in Cq. Suppose that there 
exists a nonzero function V € Cc,ia2 satisfying equation (4.24), we then claim that 
the ineqaulity |r'y(^)| < (C/*)'(^) is consequently true for all r e R and aU ^ G M. 
Writing 

5 {r e M| |r\/(0| < (C/*)'(0, C e M}, 
we verify the following properties: 

1. S* is non-empty, since € S*. 

2. S is closed. Let G 5, i = 1, 2, ... and ri —f r as i —t +00, then we will have 
\riV{£,)\ < iU*y{^) which implies that \rV{^)\ < (C/*)'(0, we therefore have r e S. 

3. S is open. Let r G S", we will show that there exists a, 5 > such that 
{r-S, r + S) C S. We claim that \rV{C)\ < (?7*)'(C) implies \rV{0\ < {U*y{0- In 
fact, let W{^) = (C/*)'(0 - rV{^) then M^(0 > 0, C e M and satisfies the following 
equation: 

{Wi - cw'i + AiiWi + A12W2 = 0, 
W'^ - CW'2 + A21W1 + ^22^2 = 0, 
(wi,W2)(-oo) = (wi,iy2)(+oo) = 0, 

where v4y, i,j = 1,2 are the entries of the Jacobian ^{U*). Since A12 > 
and A21 > 0, the Maximum Principle implies that W{£,) = {wi{£^),W2{£,))'^ > 
for ^ G M, unless W is identically and the Lemma is proved. We thus have 
(C/*)'(0 - rV{^) > 0, ^ G R. Similarly we can show that (C/*)'(0 + rV{^) > for 
^ G M. The claim then follows. 

We next show that the claim further implies |fV'('C)| < (t^*)'(C) long as r 
is sufficiently close to r. According to condition (4.18) and the assumption that 
V G Ccr^CTa; for any fixed f G M, there exists > sufficiently large such that 
(e'"i« + e-'^^€)[(C/*)'(0-ry(0] > for all C £ (-cx),iV]. This implies that fy(C) < 
(C/*)'(i^) also holds there. Furthermore, due to the claim, and the boundedness of 
the functions (U*)' and V, we can find ^ > such that for any f G {—S + r,S + r), 
we have {U*Y{£,) > fV{£,) on the finite interval [-N,N]. 

Now we fix f = f and show [(C/*)'(0 - rV{^)] > for ^ G [iV, +cx)). Noting 
the diagonal entries of the matrix ^{U* {+00))) are both negative, we can choose 
column vector P+ > such that (increasing N if necessary) ^(C/*(+oo)))P+ < 
for^G [iV,+oo). 

We have to consider the following two cases: 

Case A. If we already have [([/*)'_- fV{^)] > for f > A^, then the Maximum 
Principle impfies that [(f/*)'(C) — ^^iO] > on [A^, +00). Analogously, we have 
[([/*)'(^) + fV{^)] > is also true for ^ G M. Consequently, S is open. 
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Case B. If there is a point in the interval {N, +00) such that one of the com- 
ponents of vector ~ ^^ViO takes negative local minimum at this point, 
we consider function W{^) := (C/*)'(^) — fV{£^) + tP+. The asymptotic rates 
of ([/*)' and V at +00 imply that there is a sufficiently large r > such that 
W{^) = {U*yiO - rV{^) + tP+ > for ^ G (iV,+oo). We further assume that 
one of the components of W{S^), say wi for example, takes minimum at a finite 
point ^2 in (iV, +00). It is not hard to verify that there is a r = T2 such that the 
corresponding W{S^) staisfying ^1(^2) — and M^(^) > for ^ e (A^, +00). For 
such T2 on the one hand, we have 



= r2f^(C/*))F+ < 0; 



(4.26) ~ "i^^"" '"i^-au 



on the other hand at ^ = ^2: the first component on the left hand side of (|4.26p 
is larger than or equal to zero. We then have a contradiction, and consequently 
([/*)'(^) — fV{^) > for f e (A^, +00). We are again in the situation descibed 
by case A. By a similar argument, we can show that (J7*)'(C) + ^^{0 ^ for 

Ce[7V,oo). 

In summary, both case A and Case B show that for any r e {—5 + r, S + r), 
\fV{^)\ < (f/*)'(0,^ G R, i.e., S is open. 

Now the set 5 is a non-empty, open and closed subset of R, hence S* = K. How- 
ever, this is impossible by the definition of S, since is bounded. Therefore 
the equation LV = cannot have a nontrivial solution in Ca-ia-2- 

□ 

The next lemma shows that there is no eigenvalue of the operator L in C^-^^^^ 
with positive real part. 

Lemma 4.5. Let Cq be the complexified space of Co(R) and X be an eigenvalue 
of the operator C, given by (4-16), with corresponding eigenfunction U_ £ Cq , then 
ReX<0. 

Proof. Let the eigenvalue A = Ai + X2i and eigenfunction U_{£,) — U^{£,) + iU'^{^) 

for C G R, where A, G R and G Co(R). 

Consider the Cauchy problem ( |02-Batesj . [28-Xuj ): 



(4.27) Vt=CV-X^V, F(e,0) = C/i(O. 

It is easy to verify that F(^,t) = (^) cofi{X2t) - U'^ (Cj sa\{X2t) solves (|i:?7l) for 
^ g M and t>Q and is bounded. We suppose that at least one of the components 
of V assumes positive value at some ^ and t (we can consider —Vii otherwise). 
Suppose Ai > 0, then the following claim is true. 

Claim: There exists a r > such that V{£„ t) < rT (C/*)'(C) /or ^ e M and i > 0. 
(Recall the operator T is defined in (4.15).) 

In fact since the vector T[U* )'(^) +00 as ^ ^ —00 we can choose a sufficiently 
large ^0 > such that 

(4.28) V{i, t) < r{U*)'{0 for f < -fo and t > 0. 

Furthermore the positivity of T(L/*)'(^) for ^ e M implies that there is a r > such 
that 
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(4.29) V{^,t) < rT{U*)'{0 for $ e [-^0,^0] andt > 0. 
Let f — max{l, r}, we then have 

(4.30) V{tt)<fT{U-*)'{0 

for ^ < ^0 and i > 0. We next adjust f suitably such that an equahty in (|4.30|) 
holds on at least one component at a point (fi,ii) with ^1 G (—00,^0] and ti > 0. 

We proceed to show that the assumption Ai > implies that (|4.30p is also true 
for C > Co and i > 0. 

From the limits of gi,g2, the choice of ai,a2, and (4.21), we can find a vector 
P+ > and increasing if necessary such that 

(4.31) M(C)P+<0 fore > Co, 
where 

We can also choose a small e > such that 

(4.32) ^'"'(1) -e(2,gi+c*)/ + Af(e))P+ <0 for ^> Co, 

Now suppose that we can find a d > and a > such that V{Ci,ti) > 
fT(C/*)'(a). Let 

g+(0 

Since Q+(C) +00 as ^ ^ +00, there is a f > such that V{C,t) < rT{U*)'{£,) + 
f(5+(C) for all C > Co and t > Q, and for at least one index j and a C2 > Co and a 
^2 > 0, we have the equality for the j — th component: 

F,(C2,t2) = r-r(c/;)'(C2) + rg+(e2). 

Let F(C,t) = fT{U*)'{C) + rQ+(C) - V"(C,t), then Yj has the following properties: 
'^5(6,^2) = 0, Y,(Co,i) > 0, Yj{i,t) > for all C > Co, * > and it then follows 
that Yj^tiC2,t2) = 0, Y,,^(C2,i2) = and Yj,^^{C2,t2) > and that Y,(C,i) satisfies 

= -{CV-X,V), 

(4.33) > {-£V + XiV + CfT{U*y + C?Q+-Xi{fT{U*y + rQ+))j 
= (CY-XiY), 

= Yj^^^ -{2gl + c* )r,- ^ + Ah, {U*)Yi + Ah, {U*)Y2 - AiF, . 

Note that in the third hue above we have CfT ([/*)' = 0, CfQ+ < by (4.32), Ai > 
0, fT{U*y > and fQ+ > 0. In the last hne of (4.33), Mij denotes the ij-th entry 
of the matrix M in (4.31). However, at (C2,i2) the left hand side of l|4.33p is equal 
to zero, while the right hand side is greater than or equal to zero because the off 
diagonal entries of M is nonnegative. We have a contradiction and consequently 
the claim is proved. 

We have thus established the fact that the set 



S = {r> 0|y(C, t) < rT{U*y{0 for C e K} 
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is non-empty. Let ro denotes the greatest lower bound of S. In what foUows we 
will show ro = 0. Suppose ro > 0, we have 

(4.34) V{tt)<roT{U*Y{Otoi-^eR. 

We first assume that an equality occurs at a point in (j4.34p at the i — th component 
at a point t) with ^ G M and t > 0. Let 

Xi^,t) ■.= roT{U*y{0-V{tt). 
From (4.27), we obtain the foUwing inequality 

X, t > [CX - XiX)i 

= X,,^^-i2g, + c*)X,,^ + Mu{U*)X,+M^,{U*)X2-X,X, 
> X,,^^ - (25, + c*)X,,5 + Mu{U*)X, - XiX, 

The last inequality is true because Mij{U*) > ii i ^ j. From the positivity 
theorem for parabolic equations we deduce that Xi{£^,t) > for ^ e M and t > i (cf 
p. 14 in [16-Leung] ). However by the t— periodicity of V, we have that Xi{^,t) > 
for all ^ S M and t > 0. Contradiction with the existence of ^. This shows that 

(4.35) y(C, t) < roT{U*y{^) for ^ G K, t > 0. 

Again since T([/*)'(^) +oo monotonically as ^ ^ — oo, there exist a suffi- 
ciently small (5i > and a large M > such that 

(4.36) V{^,t) < {ro-di)T{U*y(C) for ^ < -M and i > 0. 

The positivity of T([/*)'(^) implies that we can extend inequality (|4.36p to the 
interval (— oo, M) with an even smaller S > 0. We have 

V{^,t) < {ro - 6)T{U*y{^) forC G (-oo,M) and t > 0. 

We are then in the same situation as in the proof of the claim at the beginning 
of this lemma. Using similar arguments as in the proof of the claim, we extend the 
above inequality to: 

V{^,t) > {ro - S)T{U*y{0 for ^ e R, i > 0. 
It then follows that ro — S ^ S. Contradiction with the definition of tq. Hence 
ro = 0. However, this contradicts the assumption that at least one component of 
V assume positive value. Thus we must have Ai < 0. This concludes the proof of 
the lemma. □ 

Theorem 4.6. Assume [HI] to [H3J and that ai and 02 satisfy (|4.18p . the operator 
L in Cai,a2 has a dense domain of definition. For any complext number with Re X > 
large enough, (A — L)^^ exists and is defined on all of Cai,a2i o.'^^d satisfies the 
following estimate 

(4-37) l|(A-i)-ic ^T^' 

where c > is a constant. 
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Proof. The proof follows the same idea as in |24-Volpert] but with resolvent esti- 
mates in Cq.t replaced by in the space Ccri,(T2- We skip the proof. □ 



Theorem 4.7. Under the hypotheses of Theorem \4.7\ the operator L generates an 
analytical semigroup in 0^1.^2 j where <7i and satisfy l|4.18p . 



Proof. The conclusion follows from Theorem (|4.7p and Hille-Yoshida Theorem. □ 

We now state the stability theorem, 

Theorem 4.8. Assume hypotheses [HI] to [H3], and o'i,cr2 satisfy (4-18). The 
traveling wave solution U* of p.6p - p.7p . with wave speed c* > 2-\/a , is asymptot- 
ically stable according to norm || • || || • ||c„j • That is, there exists e > such 
that if the initial condition ?7(^,0) = U{£,) G C with {U{£,) — U*{£,)) G Co-i.o-a o,n-d 
\\U — U*\ \ < e, then the solution U{^,t) exists uniquely for all t > and 

(4.38) \\U{C,t) ~ U*{x^ct)\\ < Me-^\ 

Here, the constants M > 0, 6 > are independent oft and U. 

Proof. The stability of U* leads to the consideration of the stability of the trivial 
solution for system (4.2), and the analysis of the spectrum of the operator L in 
(4.3). Corollary 14.31 Lemma 4.4 and Lemma [4.51 show that the spectrum of the 
operator L in the space Cai,a2 is contained in a closed angular region in the left 
open complex plane. Thus, following the methods in Theorem 2.1 on p. 227 in 



24-Volpert , we obtain the conclusion of this theorem. □ 
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